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SOME RAYLEIGH-TAYLOR STABILIZING MECHANISMS IN 
THIN LASER TARGETS 


I. INTRODUCTION 

The Rayleigh-Taylor^ instability arises in the acceleration of 
a fluid by one of lower density. This instability represents a 
potential obstacle to achieving laser- induced fusion in that it might 
destroy the symmetry of an imploding fuel pellet. 

There are several surfaces in a laser-driven plasma at which 
the growth of Raylelgh-Taylor instability can be large. One such 
surface is the ablation layer separating the cold, high density fluid 
from the hot, low density material ablated from the surface-. In an 
Inhomogeneous target an instability occurs at any fluid interface 
where the acceleration is directed towards the denser material. 

Several stabilizing mechanisms for the Rayleigh-Taylor 
instability are well known. In this paper we consider several 
mechanisms that tend to reduce Che instability in Chin, inhomogeneous 
laser targets. Some of the mechanisms we discuss have been treated 
in Che literature. Our reason for reviewing them here is Che 
increased relevance of chin-shell pellet designs In laser-fusion 
schemes. 

Inhomogeneous pellets have been considered for use in laser-fusion 
schemes for some time. However, recent claims from Che Lebedev 
Physical Institute^ suggest that certain inhomogeneous hollow pellets 
with very thin shells having a ratio of radius R to thickness AR of 
as much as 100 can be imploded to yield an energy gain of as much as 
10^. The number of e-foldlngs during time t of a perturbation of 
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wavenumber k expected In a homogeneous, incompressible, inviscid 

Viatfina nn a Moalnot^ la A 


material having acceleration g against vacuum is (kgt in the 
standard, linear Rayleigh-Taylor problem.^ The mode most likely to 
break up an imploding shell has wavenumber k ~ 1/AR. Then if the 


time of implosion is given by R ~ gt^/2, the number of e-foldings 

-te r>f Th-ta o-iirla 


is of order (2R/AR) . This crude analysis indicates that shells 
having an aspect ratio R/AR < 100 should be highly susceptible to 
Rayleigh-Taylor instability. Claims to the contrary warrant further 
consideration. 


The stabilizing mechanisms discussed here reduce the growth rate 
of the Rayleigh-Taylor instability in thin shells but probably cannot 
by themselves account for stability of shells having aspect ratio of 
100 . 


II. Finite Layer Effects 


In seeking effects that tend to stabilize imploding thin shells 
against the Rayleigh-Taylor instability, we consider first the 
simplest case, and find a rather large effect. We examine the effect 
of the finite widths of the fluids on the growth rate of perturbations 
at the fluid Interface. The simplest case is represented by a model 
consisting of two slabs of inviscid, incompressible fluid of finite 
widths, d^, d^, and uniform densities, and p^. 


The linearized two-dimensional fluid equations for inviscid. 


incompressible fluids at rest are' 


— Su + ik6v = 0 
dx 


(divergence-free velocity) 


v«p + “ 0 


(conservation of mass) 


VP«u = - ~ «P + gdp 


(conservation of x-momentum) 


vp6v = - lk6P 


(conservation of y-momentum) 
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Here x6u + y6v is the velocity perturbation, 6P the pressure 
perturbation and 5p the density perturbation of a fluid element of 
density p, subjected to an effective gravity gx. The perturbations 
have been decomposed by a Fourier transform into normal modes, whose 
dependence on y and t is given by 

exp (vt + iky), 

where v is the constant growth rate, and k is the constant wavenumber 
of a Rayleigh-Taylor mode. 

The equations (1) reduce to 

-h = “7 * < 2 ) 

At a fluid interface the normal component of velocity is continuous, 
as is the quantity found by integrating (2) across the interface; 
the jump conditions at a fluid interface are therefore 

A[6u] = 0 (3a) 



If the two incompressible slabs of uniform density are confined 
between rigid planes at x = -d^ and x = d^, then the boundary 
conditions there are 6u = 0. These boundary conditions, together with 
the jump conditions (3) at the interface and the differential 
equation (2), which reduces to d^6u/dx^-k^Su - 0, imply a growth 
rate given by^ 


£ Pr Pa 

gk ~ Q ctnh kd + p ctnh kd (4-) 

1 1 '^2 2 
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If instead of being confined between rigid planes, the surfaces 
are free^, then the boundary condition d6P/dt = C pertains. This 
condition leads to boundary conditions identical to the jump condition 
(5b), namely 

d 

— 6u ^ 5u = 0 (free boundary) (5) 

dx Y 

With the free boundary conditions, the growth rate of Rayleigh- 
Taylor modes is 



gk ~ 0 ctnh kd + o ctnh kd 

1 2 2 1 


(5) 


The growth rates in (^) and (6) are both less than the corresponding 
rates for semi -infinite fluids, and can be significantly less for 
perturbation wavelengths long compared to the widths of the slabs. 

In the limit of small kd, the growth rate given by (6) approaches 

zero as V -• k [ gd^d^(p^-p^)/(p^d^ + The growth rate of (4) 

also scales as k in the long-wavelength limit. If for example, 

kd , and p /p are .1, then y decreases to of the thick shell 
2 2 1 
value. 

The growth rates found for this planar geometry may be applied 
to an imploding spherical shell as well by the presciption 
k -• (-L + ^)/R, where R is the radius and t is the degree of a spherical 
harmonic In a decomposition of the Raylelgh-Taylor perturbation. The 
correspondence Is valid until spherical convergence effects near the 
origin can no longer be neglected, that Is, until (dR/dt)® « R d^R/dt^ 
Is no longer satisfied.'^ 
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III. Thin, High-Density Barrier 


Some pellet designs call for a thin layer of high-density material 
to be sandwiched between the deuterium-tritium fuel shells on the 
inside and a moderate-density ablator on the outside. In this section, 
as in the preceding, we use the simplest model that provides a 
qualitative understanding of the effect of interest. In this case, 
we seek the effect of an incompressible layer of uniform density 
and width fi sandwiched between semi-infinite incompressible fluids 
of uniform densities and p^ on either side. With the boundary 
conditions of vanishing perturbation at infinity, and the jump 
conditions (3) applied at the two interfaces, the growth rate is given 
implicitly by 


(r + 





in which s (v^-k.g)/(v^ + kg). 


(7) 


We wish to consider the limit of (7), in which the intermediate 
layer is thin, 4 « 1/k, and of high density, p^ » p^ and p^ » p^, 
but in which the total mass per area p 4 remains finite. To lowest 
order in k6 the growth rate of perturbations in this model is 


^ P, - P 


Cl + 


^k4p^p^p. 


gk 'P ,+ P 




P,)' 


kfi « 1, 0^ 4 P^ 


( 8 ) 


This result indicates that an intermediate fluid layer has 
negligible effect on the stability of Raylelgh-Taylor modes of wave- 
length much longer than the layer is wide. Since the hlgh-denslty 
layers can be much thinner than the shell as a whole, we conclude that 
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for important Rayleigh-Taylor modes the high-density thin shells have 
negligible effect on stability. Moreover, if the ablator has greater 
density than the fuel, (p^ > p^), then there is no interface 
instability. 

IV. The Perturbation in the Ablaticn 

If the adverse density gradient giving rise to the Rayleigh- 
Taylor instability does not occur at a contact discontinuity, but 
occurs rather more gradually over some finite distance, as at the 
ablation front, then an important feature in addition to the growth 
rate of the instability is the spatial location of the perturbation. 

In a laser-induced implosion of a pellet, it is obviously more 
desirable to have the instability peak farther from the ablation 
surface where convection,^ f irepolishing^ compressibility, viscosity,^ 
and heat conduction® exert greater stabilizing effects, than closer 
to the pellet surface, where it is capable of causing greater damage. 
The nonlinear growth of a perturbation is more likely to break an 
imploding shell if the peak growth in the ablation occurs closer to 
the shell surface. 

In this section we consider the simplest model of the ablation 
of an implosion in order to reach general conclusions regarding the 
location of the instability. The problem has been treated elsewhere,'^ 
but we reexamine it here for its significance to the question of 
stability of laser-fusion pellets. We will identify three distinct 
stabilization mechanisms associated with the density gradient of the 
ablation. 

Our model consists of an incompressible, inviscld fluid with an 
exponential density gradient opposed to the effective gravity, so 
that the configuration is unstable. The fluid is contained between 
a free surface at x s 0, where the density is a maximum, and a rigid 
surface at x = d. The configuration is shown in Fig, 1. There is no 
convection or motion of the fluid in the unperturbed state. With 
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a density profile oaexp (-x/h) , where h is a constant scale height, 
the differential equation (2) becomes 

>''<1 - h ^>] “ e ) 

With the free-surface boundary condition (5) applied at x =0, 
and the rigid-surface boundary condition, ^u = 0, applied at x = d, 
the growth rate is found implicitly as the positive roots of the 
transcendental equation 

tanh (n^kd) = o^khPk^h'^Ccp^ (10) 

in which co s ri _ (g/h''^) + (2kh) A discrete spectrum of 

n n " 

growth rates corresponds to the real and imaginary roots 5p^ of 
Eq. (10). This much of the problem has been treated in Ref. 4. 

This problem is just a standard Sturm-Liouville eigenvalue 
problem typically encountered in quantum mechanics. Here gk/V^ plays 
the role of the eigenvalue. The eigenfunction corresponding to the 
eigenvalue gk/v^ is 

sinh CP k(d-x) 

^u = — r-c — exp (x/2h) (11) 

n sinh CP kd r \ / 

n 

If cp^ is imaginary, as it is for all k a l/2h and all excited-state 
eigenfunctions, then (11) can be written 

sin u, k(d-x) 

6u = exp (x/2h) 

” sin u kd 
n 

in which u = icp . 

n n 

Since the growth rate is given by 
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_n. _ 4kh 

gk 1 + (2kh)‘^(l-<p“) , 

n 

the ground-state i.^genvalue (highest growth rate) corresponds to the < 

highest allowed val e of Maximum growth rates are plotted versus 

slab width for several Rayleigh-Taylor modes in Fig. 2. If the 
wavenumber k s l/2h, then the maximum growth rate is less than (kg)^ 
for all slab widths, and is given by 


max 


4kh 

(2Kh)^ + 1 




( 12 ) 


Ground-state eigenfunctions, *u, for several Rayleigh-Taylor 
modes on a given density configuration are shown in Fig. 1. The 
ground state eigenfunctions have no nodes, but excited state wave- 
functions 'SC I Hate sinusoidally within exponentially growing envelopes. 
That the eigenfunctions peak away from either boundary is a consequence 
of the boundary conditions, and will be the case for almost any 
boundary conditions. In the limit that the slab width becomes very 
large, the solutions become insensitive to boundary conditions, and 
the eigenfunctions peak towards the low density end of the 
configuration. 

Three distinct stabilization mechanisms can be identified in 
Figs. 1 and 2. first, we have seen in (12) that perturbation 
wavelengths less than the scale height have reduced growth rates 
because they "sense" only a limited region of adverse density 
gradient. This effect is illustrated in the lowest curve in Fig. 2 
for which kh = 1. Second, a pronounced stabilizing effect for all 
wavelengths is also seen in Fig. 2. The growth rate is reduced 
significantly from its maximum value if the slab width is less than 
a few scale heights. This is a thin-layer effect like that discussed 
in Sections II and III, 
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Finally, for any continuous density profile, the perturbation 
cannot peak at the maximum in density, as can be seen from Eq. (2), 
and as is illustrated in Fig. 1. In fact, it can only peak 
where the adverse logarithmic density gradient is sufficiently 
steep that -d^,np/dx > Y^/g, since that is the condition for 
d^6u/dx^ < C and d6u/dx = 0. This suggests not only that shorter 
wavelengths are stabilized in an actual implosion, but that the 
perturbation is greatest somewhere outside the surface of maximum 
density. Irrespective of boundary conditions, the perturbation is 
most likely strongest near the maximum of - d-t-np/dx. The farther 
from the pellet surface the perturbation peaks, the greater will be 
the influence of stabilization mechanisms such as viscosity and 
convection. 

We note, however, that a more complete analysis is necessary to 
determine the full extent of the effects of such stabilizing mechanisms. 
Including these effects will alter not only the growth rates but the 
spatial character of the perturbations as well. Particularly in its 
nonlinear growth, a perturbation alters its medium and the boundaries 
of its medium, thereby altering the modes of its growth. Nonlinear 
growth may cause bubbles in the density to form and break off.® "This 
process might eliminate certain modes and create new ones. 

V. Conclusions 

We have discussed in this paper several effects that tend to 
reduce the severity of the Rayleigh-Taylor instability for thin, 
inhomogeneous, incompressible, and inviscid laser targets. Some 
of these mechanisms have been treated elsewhere,*^ and are reviewed 
here for their relevance to the possibility of imploding thin shells 
stably. The effects not^.d here are probably insufficient by 
themselves to permit stable implosion of very high aspect ratio shells. 

The authors wish to express their gratitude to Dr. D. L. Book 
for helpful conversations. 
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X- 


Fig. 1 — Straight line represents exponential density profile in arbitrary units with 
scale length h and width d = 6h. Curves are longitudinal velocity perturbation 
eigenmodes 5 Uq in arbitrary units for the three wavenumbers indicated, assuming 
free surface (constant pressure) at x = 0, and rigid surface at x = d. Eigenmodes 
are ground-state, having maximum growth rates. 
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Fig. 2 — Square of Rayleigh-Taylor growth rate in units kg versus slab thickness 
in units 1/k for the density configuration of Fig. 1 and the wavenumbers indicated. 
Growtl. rates shown are the maximum for each wavenumber, corresponding to the 
ground-state eigenmodes. 
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